extern int dscal_(integer *n, doublereal *da, doublereal *dx, integer *incx);
